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1 Introduction 

A curious connection between certain integrable quantum field theories and the theory 
of the Schrodinger equation has been the subject of some recent work P}-ffl|. In this paper 
we extend these results by establishing a link between functional relations for ^-related 
Bethe ansatz systems (see for example P,|[) and third-order differential equations. Most 
of our analysis concerns a certain specialisation of the model, a particularly symmetric 
case that can also be related to the dilute A-model of j?]]. 

In the cases studied in |j}-|[], the most general differential equation was a radial 
Schrodinger problem with 'angular momentum' I and homogeneous potential x , ini- 
tially defined on the positive real axis x S (0, oo) : 

* + X 2M + l l+R) ^ x E ) = E1>(x, E) . (1.1) 

dx A x z / 

The relevant integrable quantum field theories were the massless twisted sine-Gordon 
models or, equivalently, the twisted XXZ/6-vertex models in their thermodynamic lim- 
its, and their reductions. It is worth noting that these models are all related to the Lie 
algebra A\. Spectral functions associated with ( [Cj] ) satisfy functional relations [||, [|, 
and these were mapped into functional equations appearing in the context of integrable 
quantum field theory in |l|-Q. We will follow a similar strategy here, taking a simple 
third-order ordinary differential equation as our starting-point and showing that the 
Stokes multipliers and certain spectral functions for this equation together satisfy rela- 
tions which are essentially the analogues, for the Bethe ansatz systems treated in |||7|], 
of the T-Q systems which arise in the context of the integrable quantum field theories 



related to A% [1C, 11]. This is the subject of §2, while in §3 we borrow some other ideas 
from integrable quantum field theory in order to derive a nonlinear integral equation 
for the spectral functions, an equation which is put to the test in a simple example 
in §4. Duality properties are discussed in §5, allowing us to find the equivalent of the 
angular-momentum term in ( |1.1| ) for the third-order equation. Connections with various 
perturbed conformal field theories are discussed and tested in §6. Finally §7 discusses 
the most general ^-related BA equations that arise in this context, and §8 contains 
our conclusions. 

2 The differential equation 

We begin with the following third-order ordinary differential equation: 

y"'{x,E)+P(x,E)y(x,E)=0, (2.1) 

and initially restrict ourselves to purely homogeneous 'potentials' x 3M , giving P(x,E) 
the form 

P(x, E) = x 3M - E . (2.2) 



These are the simplest higher-order generalisations of the I = cases of ( |1.1| ), and so 
we expect that some of the properties of that equation, used in the analysis of Q, will 
be preserved. In particular, motivated by the results of M for second-order equations, 



we suppose that (|2.1| ) has a solution y = y(x, E) such that: 

(i) y is an entire function of (x,E) though, due to the branch point in the potential 

at x = 0, x must in general be considered to live on a suitable cover of the punctured 

complex plane; 

(ii) y ,y' = dy/dx and y" = d?y/dx 2 admit, for M > 1/2, the asymptotic representations 

,, I T A/+1 I M+l ,. I A/+1 

y~x" A/ e M+i x , y'~- e m+i x ; y " „ x M e M+i x ; (2.3) 

as x tends to infinity in the sector 

47T 

I argx I < . (2.4) 

1 s ' 3M + 3 y J 

Furthermore, these asymptotics, or even just the asymptotic of y{x, E) with x remaining 
on the positive real axis, characterise y uniquely. 

For M < 1/2, the story is complicated by the appearance of extra terms in the 
asymptotic fl2.3p. The behaviour of the solution which decays as x — ► +00 can be more 
generally found from the formula 

rx 

y(x,E) ~ P(x,E)- 1/3 exp(- / P(t,E) 1/3 dt) , (2.5) 

Jx 

with the constant xo being related to the normalisation of the solution. (This is the 
analogue of an approximate WKB solution of a Schrodinger equation.) Since to take 
M < 1/2 would bring other technical problems into the treatment to be given below, 
from now on, unless otherwise stated, we shall restrict ourselves to M > 1/2. This 
range is the analogue of the 'semiclassical domain' of ||ll| (see |l], ||, || for a discussion 
in the context of differential equations). 

Given y(x, E), bases of solutions to the third-order equation can be constructed just 
as in the second-order case. For general values of k, define 



with 



Then y/% solves 



y k (x,E)=^y(^x,u;^ M "E), (2.6) 



, 2iri , 

u = exp — . (2.7) 

' l 3M + 3 ' 



y%(x, E) + e' 2km P(x, E)y k (x, E) = , (2.8) 

and so when k is an integer it provides a (possibly new) solution to the original prob- 
lem ( |2.1| ). However, since we will shortly need to consider fractional values, we will 
leave k arbitrary for now. It is convenient to define sectors S k as 

2&7T 



arg x 

s 3M+3 



< 3Frt • (2 - 9 > 



On the cover of the punctured complex plane on which x is defined, the sector S k 
abuts the sectors S k -\ and Sj. + i, and the sector (|2,4[) is 5-3/2 U 5- 1/2 U S\/2 U ^3/2- 
The pattern of dominance and subdominance of solutions is more involved than in 
the second-order case, since there are now three different behaviours for solutions at 
large |x|. In addition to a solution with leading behaviour x~ exp(— x +1 /(M+1)) 
as \x\ — ► +00, there are also solutions which behave as x~ M ex.p(e zizm ' 3 x M+l /(M+l)). 
(This is simply a consequence of the fact that the three third roots of —1 are —1, e 7 ™' 3 
and e -7 ™' 3 .) Depending on the sector, either one or two of these solutions tend to zero 
at large \x\. We call 'subdominant' the solution which tends to zero fastest in a given 
sector; then, up to a scalar multiple, y k is characterised as the unique solution to (2 
subdominant inside S k . 

The asymptotic (|2.3|) and the definition (|2.6[) together imply 



Uk ~ v 



(M+l)fe -M -]w^- (M+1) ^ M+1 



y k e M+i 



UJ -(M+l)k x M+l 



11 -<M+l)k M - iifi-i 



I ~{M+l)k M+l 



for \x\ — > 00 with 



x e 5 fc _ 3/ 2 U S, 



5 fc-l/2 



U 5fc+l/2 U Sfc+3/2 



(2.10) 



(2.11) 



Comparing y k , yk+i and y/c+2 in the region 5/ c+1 / 2 U5fc + 3/2, where the asymptotics of all 
three are given by ( 2.10J ), establishes their linear independence. The set {y k , 2/fc+i, yfc+2} 
therefore forms a basis of solutions to ( |2.8[) (and, for k integer, to ( p.l[) ). Alternatively, 
we can examine 

W kl ,k2,k 3 = W[y kl ,y k2 ,y k3 ] , (2.12) 

where the generalised Wronskian W[f,g, h] is defined to be 

7 /' /"' 

9 

h 



Det 



9 
h' 



9 
h" 



(2.13) 



It is a standard result (see, for example, [12] ) that, for /, g and h solving (2^), W [/, g, h] 
is independent of x, and that /, g and h are linearly independent if and only if W[f, g, h] 
is nonzero. For (k\, k%, ^3) = ( — 1, 0, 1), the asymptotic ( P-IOQ , used in iS-i^USi^, shows 
that 



W- 



1,0,1 



2M 



^ sm ( 3 M+3 7r ) sin ( 337+3^ 



It is also the case that 



Wk-L+afo+ato+aiE) = W^foM^ C '^O 



(2.14) 



(2.15) 



so Wk,k+i,k+2 is nonzero for all k, thus confirming the independence of {y^, yk+i, Vk+2}- 
We now aim to generalise the analysis of [§| to this situation, guided in part by the 
treatment of ^-related BA systems provided by ||. Since yi,y2>2/3 form a basis, we 
can write 

y - S (1) (£) yi + 5( 2 ) (E)y 2 - y 3 = (2.16) 



with 

s (i) {E) = Wo^ ^ s (2) {E) = WifiA_ (217) 

Wl,2,3 ^1,2,3 



The coefficient of 1/3 in ( |2.16| ) is —1 by ( 2.15 ); S^- 1 ' and 5"' 2 ' are Stokes multipliers for 



( |2.1| ), and are analytic functions of E. Notice the formal similarity between this equation 
and eq. (15) of ||. 

Now suppose that k\ and k 2 differ by an integer. Then y kl and y k2 both solve (2.8) 



(with e 2km = e 2kl7Tl = e 2fc27n ), and it can be checked by direct substitution that the 
function 

z klk2 (x,E) = y kl y' k2 - y' kl yk 2 (2-18) 

solves 

4"fc 2 (*, #) " e" 2fcm P(x, P)z fclfc2 (x, E)=0. (2.19) 

This is just the equation adjoint to ( |2.8D ; the observation that the Wronskian of two 
solutions of a third-order ordinary differential equation satisfies the adjoint equation 
dates back at least to Birkhoff jy! . Observe also that if k\ and k 2 are shifted by a 
half-integer, then a solution of the original equation fl2.8Q results: 

z>" 1 i^fij+e-^Pf^K , 1, , i(x,£) =0 . (2.20) 

fcl + 2 fc 2+2 K1 + 2K2+2 

For |fci — &2| < 3, the regions ( |2.1l| ) for k = k\ and k = k2 have a nonempty overlap, 
and an asymptotic for Zfe 1 fe 2 is easily obtained from ( 2.10| ). In particular, for k = 1,2,3 
we have 

Z-k/2,k/2(x,E) ~ 2isin(7rA;/3)x" M e~ 2cos(W3) A7+T a ' M+1 , s -» +00 . (2.21) 



For k = 1, 2_i/2,i/2 solves fl2.l|) , and now from ( 2.21 ) we see that it shares (up to a 



proportionality factor) the asymptotic (pO|). By uniqueness, we deduce 

z. 1/2 ,i/2(x,E)=iV3y(x,E). (2.22) 

Unfortunately, this argument is not so effective for the other cases. At k = 2, the formula 



( |2.21 ) shows only that 2-1,1 is not subdominant on the real axis, and this information 
is not enough to pin the function down. For k = 3, sin(7r£;/3) = and all that can be 
deduced is that the leading asymptotic of £-3/2,3/2 i s subleading to the term over which 
we have control. 

The next step is to manipulate Q2.16 ) in order to eliminate either S^ or S^ 2 ' . We 



have 

y , 1 y -S^(E)y' 1 y 1 + S^(E)y' 1 y 2 -y' 1 y 3 = 0; (2.23) 

y 1 y' -SW(E)y 1 y' 1 + S^(E)y 1 y' 2 -y 1 y' 3 = 0, (2.24) 

and, subtracting, 

S^(E)z 12 = z 01 + z 13 . (2.25) 



For the reasons just explained, functions Zk x k 2 with \ki — k 2 \ = 1 are the most easily 
handled, so we use the identity y 2 z±3 = yi-^23 + #122/3 to rewrite ( |2,25 ) as 



S^\E)y 2 zi2 = y2Zoi + z 12 y3 + ^232/1 • (2.26) 

Likewise, 

5 (1) (E)yiz 12 = z u y + yi^23 + ^012/2 • (2.27) 

Now the result ( J2T22] ) can be combined with shifts in E to w 15M / 4 £ and u? XM I^E 
respectively to rewrite both ( 2.26 ) and ( |2.27| ) as 



T{E)y_ l/A y l/A = y_i/42/5/4 + 1/-3/42/3/4 + 2/-5/42/1/4 , (2-28) 

where 

T(E) = S« (cu 15M / 4 £) = S< 2 > (w 21M / 4 £) . (2.29) 

As a byproduct, this has established that the two Stokes multipliers S^- 1 ' and S^ 2 ' are 
related by an analytic continuation in E. 

Finally, taking Q2.28 ) at x = yields a functional relation involving E alone. To 



absorb various phases, it is convenient to set 



Q+(E) = E-*My(0, E) , Q+(E) = Q + (uT 3Mfe £) . (2.30) 
Then the relation is 

T Q-l/4:Ql/4 = Q-I/4Q5/4, + Q-3/4^3/4 + Q-5/4^1/4 ' (2.31) 

This is very similar to the equations related to the dilute A model studied in [0,[l4j] . An 

equation involving y"(0,E) can also be derived. First, differentiate (|2.25| ) twice with 
respect to x: 

S^(E)z'{ 2 = 4 1 + z'{,. (2.32) 



Using the fact that y\, y 2 and 773 all solve ( |2.1| ), we have y 2 z'{^ = y'{z 2 ^ + Z\ 2 yi,i an< ^ so 
the previous steps can be repeated to find 

T(E)y': i/4 y'l /4 = y'U/^/A + ^-3/4^3/4 + ^U/^i/A ■ (2-33) 

Again set x = 0, and define 

Q-(E) = ±Etmy"(0,E) , Q~ k {E) = Q-(io- 3Mk E) (2.34) 

(the factor * is included for later convenience). Then 

^Q-l/4 ( 3l/4 = Q-I/4Q5/4 + Q-3/4^3/4 + Q-5/4^1/4 ■ (2.35) 

There is no simple relation involving y'(0, E) alone, but from ( [2.18D and ( |2.22| ) one can 
deduce 

iV3y' = y-i/ 2 y'l/ 2 - y'U/2Vl/2 , (2-36) 

which allows y'(0,E) to be recovered once y(0,E) and y"(0,E) are known. 



3 The non-linear integral equation 

The functions Q^(E) are not single- valued, and to derive an integral equation it is more 
convenient to work with the functions y(0,E) and y"(0,E) directly. Set 

D+(E)=y(0,E) , D-(E) = ±y"(0,E) (3.1) 

(so that D ± (E) = E ± ^iQ ± (E) and Q^(E) = u ±k E T ^T D ± (uj- 3Mk E) ). These are 
entire functions of E and can be interpreted as spectral determinants for the third- 



order equation (2.1), since their zeroes coincide with the values of E for which the 
solution y, decaying at x — ► +oo for all values of E, in addition either vanishes at x = 
(for the zeroes of D + ), or has a vanishing second derivative at x = (for the zeroes 
of D~). (See, for example, ref. Q for a more detailed discussion of this point in the 
context of second-order equations.) For M > 1/2, the functions D (E) have large- 1 .27 1 
asymptotics 

lnD ± (E) ~ a (-E) ft \E\ -> oo , |arg(-£)| < vr (3.2) 

where fi = (M + 1)/3M, a = k(3M, 3), and 

The growth of In D(E) is no larger on the positive real E'-axis than elsewhere, so the 
orders of D + and D~ as functions of E are both equal to /i, and are less than 1 for 
M > 1/2. Invoking the Hadamard factorisation theorem, we can write 



^■^n('-|). 



(3.4) 



The precise values of the constants D (0) are irrelevant for the treatment below, but 
some knowledge of the positions of the zeroes {E k } will be crucial. We conjecture that, 
for all M > 0, all of the zeroes of D (E) lie on the positive real E-axis. Some numerical 
evidence in favour of this claim will be presented below. 

The generalised T-Q relations (^3l|) , fl2~35l ) taken at either E e {co 3M/4 E±} or 
E e {uj- ui /*E±} imply 

" " W n+, W/9rn±N ' ( 3 ' 5 ) 



D±(w- 3M /2£±) " D±(w 3M / 2 ^) 

an equation that can be written in a Bethe-ansatz form as 

°° P ± , -3Af p± °° p± , .-3M/2 p± 

11 B±_ W 3M £ ± " W 11 £±_ w ^/ 2 ^ • ^ 

This equation is at least not inconsistent with the conjectured reality of the E n 's, since 
both sides then reduce to pure phases. There are certainly other, complex, solutions to 



( p.6p , so the reality property should be seen as a way of selecting the particular solution 
relevant to our differential equation, analogous to the selection of the ground state in 
an integrable model. 



A non-linear integral equation, similar to those described in [11, 15, In], can now be 
obtained for the quantity 



± ±1 D± (u- 3M E) D ± (u 3M / 2 E) 



d ^ E )- W D ±( U 3M E) D ±( U -3M/2 E) ■ ( 3 - 7 ) 



We shall follow a path that completely parallels the treatment given in pi] . By 
d ± (E) = —1 at the points {E k }. (The value —1 might also occur at other points; we 
supplement our previous conjecture with the assumption that none of these points lie 



on the positive real axis.) The product representation (3.4) implies 



2 
In d ± (E) = ±m^— - + Y, HE /E n ) (3.8) 

n=l 

where 

F(E) - ln (1 - ^- 3M ) (1 ~ S" 3M/2 ) (39) 

' ^ j " (1 - ^3M) (1 _ ^-SM/2) • ^-9) 

The sum over the E n in ([T^) can be written as a contour integral 

ln^(£) = ±^^— - + jf _F(£/£')^ ln(l + rf*(tf)) (3.10) 

with the contour C running from +oo to above the real axis, winding around and 
returning to +oo below the real axis. (It is at this point that the conjectures about the 
locations of the E n 's and of the other zeroes of d (E) + 1 are used.) If the new variable 
9 = filnE is introduced, the function F becomes 

with £ = 1/M. Now define 

f ± (e)=lnd ± (e 3Me / {M+ V), (3.12) 

use the property d ± {E)* = d ± {E*)~ 1 and integrate by parts to recast (|3.10|) as 

/OO 2 

d0'i?(0 - 9') ln /±(0' - iO) = ±*tt^— - 

/OO 
d0'.R(0 - 003m ln(l + Z^' - *0)) (3.13) 

-OO 



with 

ij(0) = J_5 e F(e 3M0 /( M+1 )) . (3.14) 

2-7T 

The term (1 — R) * In/ (0) on the LHS of ( |3,13 ) is easily inverted using Fourier trans- 
forms. Using 

smh(h6 + iirr) 2/isin(2r7r) 

sinh(/i# — i7rr) cosh(2/i(9) — cos(2t7t) 

dfl ifcfl 2fe S in(2r7r) = sinh((l - 2r)fg) 

2vr cosh (2/i0 )- cos (2rvr) sinh(f^) ' 



we have 



sinh(f(l-Ofc) sinh(f^) _ 2sinh(^) cosh(f (l-2^)fc) 

sinh(f (1+0*:) + sinh(f (1+0*) ~ sinh(f(l+£)&) 

sinh(f^)cosh(fA;) 



1 fl(Ai) ~sinh(f(l+Ofc)cosh(ffc) " (3-17) 

Transforming back to space and rewriting the imaginary part in terms of values above 
and below the real axis, the functions / (0) solve 

/+(#) = ±ma-ib e e + [ <p(9-6')]n(l+e f± ^)d9' 

~ [ <p(9-6') ln(l+e- f±{e,) )d9' (3.18) 
Jc 2 

where a = 2/3, the contours C\ and Ci run from — oo to +oo, just below and just above 
the real 0-axis, 



/ e muu^ cusi H g m+ - z^; qg 1 



oo e ik8 sin h(A;|) cosh(f fc(l - 20) dfc 
cosh(f fc) sinh(/c§£) 2^ 

and the constant &o = 2 sin(7r J ^ j ^ 1 )ao has been fixed using the asymptotic behaviour (|3.2|) . 
(The corresponding zero mode can be traced to the zero in 1 — R at k = i .) The param- 
eter a in ( 3.1§| ) is analogous to the chemical potential (or twist) term in the equations 
of 0|l|g^ 

A first consistency check is immediate: in the large 9 limit of ( 3.1§| ) the driving 
term 6oe dominates, and so in this limit the functions exp(/ (9)) are —1 at the points 
9 = 9±, or 

E = E± = e e ^ = ((2n-l±l)n/b )» , (n = l,2,...). (3.20) 

The same limit can be treated directly using a WKB-like approach to the differential 



equation (|2.l| ). Start from (2J3) with x > xq and fix xq to be at the inversion point 



P(xq,E) = (xq = E 1 ' 3 ^ 1 ). Now, using analytic continuation (see for example §47 
in JT3I), the dominant part in the region x < xq is 



1 / rxo 1 \ /a/3 /"^o 1 

y(x,£) ~ |P(x,S)|"3exp( i / |P(x,S)|3dx)cos ^— / \P(x,E)\3dx - § 




Thus to have y(0,E) = 0, requires 

^3 r°(x 3M -E+)3=6o(S+r = (2n-i)7r , n = l,2,..., (3.22) 

JO 



'0 

where the formula 



1 1 sinf- + -") 

(l-x a )b= {b l-\ (a,b), (3.23) 

sm fe- 



was used. The prediction ( 3.22J ) agrees perfectly with ( 3.2C| ). In figure 1 the positions of 



the lowest zeroes of D + {E) are plotted in the range 0.1 < 2>M < 7, and compared with 
the WKB-like prediction. Evidence for the reality of the E n at 3M = 1 will be given in 
the next section; in the meantime, we note that the levels continue smoothly away from 
that point, and the eigenvalues appear to remain real in the range studied. (The figure 
can be compared with figures 1 and 2 of ||], which illustrate cases where the spectrum 
does not remain real in the full range displayed.) The kernel <p{0) given in ( [3. 19 ) 



coincides with i/2ir times the logarithmic derivative of the scalar factor in the Izergin- 
Korepin S-matrix for the a 2 model [[Iqj (cf. eq. (3.21) of [19|, though note that the 



normalisation of £ used by Smirnov in [jy| differs from ours: £[ Smirnov l = _iL^[this paper] ^ 
This is an element of the advertised link between the differential equation (|2.1|) and 



(2) (2) 

the a 2 model, the parameters being related as M = l/£ (with £ related to the a 2 



coupling 7 as £ = 7/(271"— 7) |19| ). When 3M is an integer the potential is analytic, and 
the associated scattering theory is diagonal; the same phenomenon was observed in the 
Schrodinger/sine-Gordon case in Q. The similarity between the relations ( 2.31[ ), ( [2.35 ) 



and ( p.6| ) and those arising in the dilute A model |7l[L4j has already been mentioned. 
Since the a 2 model is conjectured to be the continuum limit of the dilute A model (see 
for example p(|), the fact that elements of it emerge here is not a complete surprise. 
Nevertheless, it is an encouraging signal that we are on the right track. We will return 
to this point in §pl 

4 The linear potential 

A simple but non-trivial example occurs when M = 1/3, and is the analogue of the 'Airy 
case' of the second-order problem, discussed in p[pi|. This lies outside the M > 1/2 
zone treated so far, so we have to assume that the results obtained above continue to 
hold as the region of their initial derivation is left. The basic differential equation is 

y'"(x, E) + xy(x, E) = Ey(x, E) . (4.1) 




Figure 1 : The positions of the first eight zeroes of D + (E) , plotted on a log scale. Dotted 
lines show the WKB-like predictions, and solid lines the results from the nonlinear 
integral equation. 



Setting y(x, E) = A(x — E), this becomes 

A'"(x)+xA(x) =0 . 
This equation is solvable via a complex-Fourier transform: 



A(x) 



e -ipx+\p± dp 



(4.2) 



(4.3) 



where the integration path V is represented in figure 2. (A curious feature of this case 
is that the function T(E) is a constant, equal to 1.) Even though the problem is not 
self-adjoint, numerical evidence suggests that all the zeroes of A(x), A'(x) and A"(x) 
lie on the negative real axis (see figure 3), and so the zeroes of y(0,E) = A(—E), and 
of y'(0,E) and y"(0,E), are positive and real. In the first columns of tables 1 and 2, 
the positions of the first ten zeroes of A(— x) and A"(— x) are displayed. 

The approximate positions of the zeroes of A(x) can be found from the WKB formula 
of the last section. As a check, we rederive them here via a saddle-point treatment of 



10 




Figure 2: The integration path. 




Figure 3: A search for the zeroes of A(x) in the complex x-plane. 
The function plotted is |A(x)|/(l+|A(a;)|) , A(x) = e x A(x) . 



( |4.3[ ), The exponent of the integrand for x = —\x\ < has stationary points at 



po = i\x\3 , p± = ie 



±i2?r/3 



X 3 



Deforming the contour T so that it touches the points p±, we get 

•ZL 3 

! 3 +e4 



AHx \) „ M{e l\^e~^ + q +e ||,|4/a e W3_ jf ^ 



(4.4) 



(4.5) 



where the phases ±7r/3 are the contributions from the choice of the steepest descent 
directions, transforming the quadratic terms in the expansion near the saddle points 
into a pure Gaussian integral 



v = ''s 



3 r 00 3, ,| + 2 i 

e 2' i di = a; 3 . 



(4.6) 
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E k (Exact) 



E k (WKB) 



E k (NLIE) 



1 2.8868281617697677 2.84467 

2 5.1522519299627660 5.13866 

3 7.1303732976716514 7.12265 

4 8.9403621072563961 8.93513 

5 10.635608688272157 10.6317 

6 12.245164125544329 12.2421 

7 13.787063381394688 13.7846 

8 15.273489957153985 15.2714 

9 16.713173447810789 16.7114 



2.886828161769766 
5.152251929962763 
7.130373297671650 
8.940362107256395 
10.63560868827213 
12.24516412554437 
13.78706338139461 
15.27348995715393 
16.71317344781078 



Table 1: Zeroes ofA(-x). 



Ef. (Exact) 



E k (WKB) E k (NLIE) 



1 0.8909448213691012 0.85075 

2 3.6439018898747455 3.66124 

3 5.8171363674931786 5.82455 

4 7.7377094704565507 7.74230 

5 9.5084853382296383 9.51174 

6 11.174539990553549 11.1770 

7 12.761111499098718 12.7631 

8 14.284208133811659 14.2859 

9 15.754819677279076 15.7562 



0.890944821369104 
3.643901889874749 
5.817136367493172 
7.737709470456556 
9.508485338229636 
11.17453999055357 
12.76111149909874 
14.28420813381167 
15.75481967727909 



Table 2: Zeroes of A"(— x). 



Thus for large negative x we have 



l 3 ^ 

„4(-|x|)~2|x|-3 e 8 N3 



COS 



f\/3>|3 



(4.7) 



For x > 0, the dominant saddle point is instead at pq = —i\x\ 3 , and 



l 3 - 
A(x) ~ |x|-3 e -I N3 . 



(4.8) 



(This agrees with the general asymptotic (|2.5| ), since y(x, E) = A{x—E).) The dominant 
behaviours of A'(x) and A"(x) for x real and |x| large are 



A'{-\x\ 



-2e 



M3 



cos^ 



lVs\x\r 



A"(-\x\) ~ 2|x|lel |x|3 cos(|\/3|x|t + 



^'(|x|) 
A"(\x\ 



3 1 



' -ikl 3 



x 3e 4i 



(4.9) 
. (4.10) 



Continuing to differentiate, the general result for the approximate positions of the zeroes 
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of the m th derivative A^ m ' (x) is 

/ A 9m 4- 1 \ 3//4 

^M( X ) = : x = -\-—(2n--^-)irJ (« = 1,2,...). (4.11) 



At m = (4.11) reduces to the M = 1/3 WKB prediction (Bjgg). In tables 1 and 2 



the results from formula ( [4.11 ) are compared with the 'exact' result from a numerical 



treatment of ( |4.3| ) , and also against the results of the numerical solution of the nonlinear 
integral equation (|3.18| ). Clearly the agreement is very good. 

5 Duality and a more general chemical potential 

In this section we shall investigate the effect of a duality transformation on (2.1), anal- 



ogous to that studied in J2) for the Schrodinger equation ( |1.1|) . In the Schrodinger 
case duality maps wavefunctions for confining potentials (M > 0) to wavefunctions 
for singular potentials (—1 < M < 0), in such a way that the theories with M and 
M = —M/(M + 1) are dual, their respective spectral problems being essentially equiva- 
lent. It also changes the coefficient of the 'angular momentum' term l(l+l)/x 2 in (1.1) 



in a non-trivial way; the same phenomenon here will allow us to guess the corresponding 
term for the third-order problem (|2.1|). 



To implement the duality transformation, we begin with a Langer [22| -type variable 
transformation 

y{x) = e z u(z) , z = lnx , (5-1) 

after which ( |2.1[ ) becomes 

u'"{z) - u'{z) + ( e ( 3M + 3 )^ - Ee 3z )u(z) = . (5.2) 

The duality M — v M is now effected by interchanging the roles of the two exponentials. 
Substituting z — ► M z . 1 + In ^3 yields 

u"{z) - — L— ti'fr) + (- e 3-/(^+D _ Ee^)u(z) = , (5.3) 

where E = -(M+1) 3M /E M+1 . 

Now transforming back results in the equation 



M(M+2) 1 d 1,, , , _-*M, 

(M+l) 2 [ x 2 dx x 3 



y'" + 7;;;;;, (^iz - - 3 )y + (-^ - ^)y = . (5.4) 



As promised, the confining 'potential' x has been exchanged for a singular potential 
_ x -3M/(M+i)^ an( ^ a new ^ eriri) proportional to (x~ 2 d/dx — x~ 3 )y, has been generated. 
This motivates us to enlarge the set of differential equations under consideration to 

y'" - G ^i-- \)y + (* 3M - E )y = ° ( 5 - 5 ) 

x z dx x^ 
13 



with G a new parameter, analogous to 1(1+1) for the Schrodinger equation. Duality 



maps (5.5) to 

V" -G(lT- \)V + (-*** - % = , (5.6) 

where 

y(x,E,G) = (M+l)- 1 J E 1/3 x( 1 ^T)y((M+l)S- 1 /V /(M+1) , J B,G) (5.7) 

and 

~ M - (M+lf M ~_ G- M(M+2) 

m -"mTT ' ^-~ W+i ' G " (M+i)^ • (5 ' 8) 

Duality therefore maps the 3-parameter family ({M, E, G}) of differential equations ( |5.5| ) 
onto itself. The analysis of §2 can now be repeated for these generalised problems. It 
is convenient to write G = g(g+2) and to work mostly with g instead of G. We first 
enlarge the scope of fl2.6|) by setting 

y k (x, E, g) = uj k y(u~ k x, u~ 3Mk E, g) ; (5.9) 

then y k solves 

Vk - G{±± - ±)y k + e~^P(x, E)y k = . (5.10) 



Next we define z klk2 as in flglgp. If e ~ 2kl7Ti = e~ 2k27Ti = e ~ 2klT \ then 



4% " G(±± - ±)z klk2 - e- 2km P(x,E)z klk2 = , (5.11) 



which is the equation adjoint to ( 5,10| )n. We can recover the original problem by shifting 



A; by a half-integer, and arguing just as before we find that 

z -i/2,i/2(^, E, g) = iV3y(x, E, g) (5.12) 

and then 

T(E, o)y„i/ 4 J/i/4 = J/-1/4 2/5/4 + V-3/4 2/3/4 + V-5/4 2/1/4 , (5-13) 

where T(E,g) = S^(u l5M ^E,g) = S {2 \uj 21M '■> /4 E , g) , and S& and S^ are defined as 
in ( 2.16Q . A non-zero value of G = g(g+2) causes ( |5.10| ) to be singular at the origin, so 



simply considering ( |5.13D at x = is not an option. Instead, we expand y(x,E,g) as 

y(x,E,g) = D + (E,g) X + + D°(E,g) X o + D~(E,g) x- , (5.14) 

where {x+, Xo> X-} forms a basis of solutions defined via behaviour near the origin: 

Xi(x,E,g)~x x *+0(x x < +3 ) , * = +,0,-, (5.15) 

'Notice that the operator x~ 2 d/dx — x~ :i is by itself anti-self-adjoint; this gives some insight as to 
why it is a sensible generalisation of the x~ 2 term in the Schrodinger equation. 
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with the Aj's the roots of the indicial equation (A — 1)(A(A — 2) — g(g + 2)) = : 

A+ = -g , A = 1 , A_ = g + 2 . (5.16) 

Explicitly, the functions D + , D° and D~ are 

D + = W[y,xo,X-] D o = W[y,x-,X+] D - = W[y,x+,Xo] / g ^ 

W[x+,Xo,X-} ' ' W[x+,Xo,X-} ' "' W[x+,Xo,X-] 

with W[ X+ , xo, X-] = 2 (g+1) 3 . At 5 = they reduce to y(Q, E), y'{0, E) and y"(0, E)/2 
respectively, in agreement with the notation of earlier sections. 
Defining 

Q ± {E,g) = E^D ± (E,g) , Q± = Q ± {^ Mk E , g) , (5.18) 



the generalised T-Q relations ( 2.31| ), ( |2.35[ ) have exactly the same form as before, and 
( |3.6p becomes 

The arguments of §|| can now be repeated essentially verbatim, to discover that, so long 
as G is such that the conjectures of §3 about the zeroes of T and D remain true, 
the quantities D + and D~ for the more general differential equation ( |5.5| ) are again 
described by the nonlinear integral equation ( |3.18| ), but with chemical potential term 
now taking the value a = § (<7+l) • 

At M = 1, lo ±3M = -1, the LHS of ( g3g ) is 1 and the A 2 -related BA equation 
'collapses' onto one more closely linked with A\. This leads to a rather surprising 
equivalence between spectral problems for a Schrodinger equation with potential x 6 + 
1(1 + l)/x 2 and the third-order problem ( |5.5| ) at M = 1. For these special points, the 
quantities in this paper are related to those of H as 

dhe, 9 )\m=i oc d±( C - 3 /^,oiE 3 [4]1 , n^?)k=i = ikc-^^iEs 01 , 

(5.20) 
where £+| = | (g+1) , c = ^ M= ' pWy^ , and T\ is one of the 'fused' T-operators 
discussed in §4 of pj. 



6 A link with perturbed conformal field theory 

(2) 

We now return to the relation with the a 2 model, briefly mentioned at the end of §3. 
The analogy with results of [15,16] for the Ai-related models suggests that the quantity 



c e ff 



f ^Hl+e^^dO- [ e e ln(l+e-^))^ (6-1) 

JCx JC 2 J 



15 



should be interpreted as an effective central charge of an underlying conformal field 
theory. For general a = I (<?+l) , this would predict 

3 

cfeff = 1 - — — -a 2 . (6.2) 

M + 1 v ; 

Going further, it is natural to interpret ( ]3.18| ) and ( |6.1| ) as the ultraviolet limit of the 
following 'massive' system: 

f(9) = ma - ir sinhfl + / ip(9-9') ln(l+e-^')) d9' - I (p(9-9') ln(l+ e - f( - '^) d9' ; 

JCi Jc 2 

Ceff (r) = ^ ( f sinhfl ln(l+e /(e) ) d9 - f smh9 In(l+e _/ W)tW J . (6.3) 

(2) 

This should encode finite-size effects in the massive a 2 theory, with r = M S R, M s the 
mass of the fundamental soliton and R the circumference of the (infinite) cylinder on 
which the theory is living. (Notice that there is no need to distinguish / + from /~ any 
more, since the mapping 9 — > —9 now has the effect of negating a.) There is now a 
natural scale, which can be related to an operator <p perturbing the ultraviolet conformal 
field theory. Standard considerations [23], based on the 9^9 + i2ir 3j ^ periodicity of 



f(9), suggest that so long as a is not an integer c e s(r) will have an expansion in powers 
Q £ r 6M/(Af+i) (together with an irregular 'anti-bulk' term, irrelevant to the current 
discussion). This implies for <f) either the conformal dimensions 

Q /I f 

h 6 = hs = l (6.4) 

2M+2 v ; 

and an expansion of c e fr(r) in which only even powers of the coupling A to the operator 
(f) appear, or, alternatively, the conformal dimensions 

- , 3M . , 

h 4 , = h 4> = l- — - (6.5) 

and an expansion which sees both even and odd powers of A. 

When a is an integer, the standard considerations of p3[ may have to be modified. 
Absorbing the term ma into a shift in f(9), ( |6.3| ) becomes exactly odd under a negation 
of 9. This forces the shifted f(9) to be zero at 9 = 0, even in the far ultraviolet, and 
so long as the would-be plateau value is nonzero, it splits the plateau region into two 
pieces, each of half the previous length. As a result, the regular expansion of c e fj(r) is 
in powers of r 3M/(M+i) ^ an j no ^. r 6M/(Af+i)_ p orrnu i a ( |6.4| ) now describes the situation 
when both even and odd powers of A appear in the expansion of c c g, while for even 
powers only, the correct formula is h<h = h<f> = 1 — 3M/(4M+4). Note though that this 
plateau-splitting effect does not occur at a = 0, since for this case the plateau value of 
/ is anyway zero, and imposing /(0) = has no effect. 

As explained in |l^] (see also p5|-p7|] ) , the a 2 model, when appropriately quantum- 



reduced, should correspond to the minimal models M Pt q (with p and q coprime integers 
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andp < q) perturbed by either 12 , 02 1 or [28,2S] 0i5- With 12 the perturbing operator, 

(6.6) 



the relation with the parameter £ appearing in the kernel ( |3.19 ) is |19| 

p 2e 



q (i+e) 

Since M = l/£, the ultraviolet effective central charge for a given value of a, as predicted 
by Q, is 

Ceff = 1 - 7T« 2 • ( 6 - 7 ) 

2q 

To recover 02i perturbations one simply has to swap p and g in ( |6.6D and ( |6.7D fl9(| , 
while to find 0is, p/g should be replaced by &p/q pq| . 

For the sine-Gordon model, naturally associated with the 0i3 perturbing opera- 
tor [3C, 31], it has been observed both analytically and numerically [pT|, |15|, |16|, |32 , 33] 



that reduction is implemented at the level of finite-size effects and the non-linear inte- 
gral equation via a particular choice of the chemical potential. The similarity between 
our equations and those in []l5|,[l6| suggests that the same should be true here. To decide 



which value of a will tune (6.3) onto the ground state of the relevant perturbed minimal 
model, we demand that the ultraviolet effective central charges match up; the predicted 
dimensions of the perturbing operators, and a comparison of results at nonzero values 
of r with those obtained via the thermodynamic Bethe ansatz method, will then provide 
some nontrivial tests of the proposal. 

The effective central charge of the ground state of the theory M. pq is c e ff = 1 — 6/pq. 
Thus to have any chance of matching the vacua of the 0i2-perturbed models, we must 



set a = 2/p. The required value of /i^, namely hu = -K — 3, is then matched by (6.4) 



The value just chosen for a being a non-zero integer if and only if p = 2, (16^) will 
be the correct formula to use provided the regular parts of the ground state energies 
of the models M. vq perturbed by 0i2 expand in even powers of A for p > 3, and in 
even and odd powers for p = 2. This 'prediction' holds for all of the examples that 
we checked. We then compared numerical results for (M, a) = (4,1), (1,1), (5,3), 

(f,i) and (1,0) against the tables of ||,|| for the A^ ] (Yang-Lee), E 8 , E 7 , E 6 and 
Z?4-related TBA equations, respectively, finding excellent agreement. Swapping p and 



q in (6.7), the choice a = 2/q should capture the 02i cases. The conformal weight of 
^21 = 4^ — 3 is matched by fl6,4| ), provided the swap of p and q in ( |6^ ) is remembered. 
This time a is never an integer, and the use of ( |6.4[) is justified by the regular parts 
of the ground state energies of the 02i perturbations always being in even powers of 
the coupling A. For (M, a) = (3,3), (§>j) and (| 5 §) the results from the A± and 
A2-related TBA equations |24j and the M. 35 model |35| , p^| w ere reproduced within our 
numerical accuracy^]. Finally, replacing p/q by ^p/q in ( |6.2D (q > 2p), at a = 1/p the 
models M. pq perturbed by (/>is are recovered. This time it is ( |6.5D which predicts the 
correct value for h^, as expected given that ^15 perturbations expand in both even 

^Beware of a misprint in eq. (7) of Q: the (minus) sign before Yl ,=1 should be reversed. 
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and odd powers of A. TBA equations for a number of 0i5-perturbed models have been 
proposed in [28,35,37], but so far we have only compared ( |3.18| ) with the TBA for the 
015 perturbation of the .M37 model given in |35|]. A selection of our numerical results 
for all of the cases just mentioned is presented in table 3, together with thermodynamic 
Bethe ansatz data taken from refs. [p4|, 34 -|36[|. To facilitate the comparision, we took 



r = M\R throughout, with M\ the mass of the fundamental particle in the reduced 

(2) 
scattering theory. For A 2 and Eg this is the first breather in the unreduced theory, 

and Mi is equal to 2cos(57r/12)M s and 2 cos(37r/10)M s respectively. In all of the other 

models in the table, Mi is equal to M s . The results strongly support the claim that 

the system (6.3) encodes the ground state energies of 0i2, 02i and 0is perturbations of 



minimal models. 



Model 


(M,a) 


r 


TBA 


NLIE 


4 2) + 012 


(4, 


i) 


0.001 


0.399999735051974 


0.399999735051971 








0.002 


0.399998953903823 


0.399998953903824 


Eg + ^12 


\ 3' 


i) 


0.025 


0.499926331494289 


0.499926331494288 








0.05 


0.499705463734389 


0.499705463734387 


E7 + 012 


V2' 


h 


0.02 


0.699928050129612 


0.699928050129611 








0.04 


0.699712371203531 


0.699712371203531 


E& + 0i2 




*) 


0.025 


0.857016839032789 


0.857016839032790 








0.05 


0.856639509661813 


0.856639509661819 


Da + 0i2 


(1, 


0) 


0.01 


0.999972507850553 


0.999972507850552 








0.02 


0.999890328583463 


0.999890328583464 


Ai + 021 


a, 


h 


0.02 


0.499697279140833 


0.499697279140832 








0.04 


0.498957654198721 


0.498957654198722 


A 2 + 021 


(i 


*) 


0.001 


0.799999470103948 


0.799999470103940 








0.002 


0.799997907807646 


0.799997907807649 


M35 + 021 


(h 


I) 


0.1 


0.596517064916761 


0.596517064916762 








0.15 


0.592881408017592 


0.592881408017593 


M 37 + 015 


(h 


\) 


0.1 


0.709591770021299 


0.709591770021299 








0.15 


0.705031895238354 


0.705031895238357 



Table 3: NLIE results versus TBA data from refs. ||,|| |6| 



In previously-studied examples, equations similar in form to those for the ground 
state have been found to describe excited states (see for example |jll|, ^2 , 38|, 39| ) . We 
expect that the same will be possible here, but we will leave investigation of this point for 

future work. Finally, we remark that it would be interesting to derive a nonlinear integral 

(2) 
equation for the a 2 model directly from finite lattice BA equations. We understand 

that progress is currently being made in this direction f4T 
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7 General v4 2 -related BA equations 

In this section, we discuss the effect of adding a term proportional to x~ 3 to the differ- 
ential equation (|5.5|). The equation becomes 

y > "-G(\-f-±)y + P(x,E,L)y = , (7.1) 

x z ax x 6 

where 



„3M • ■ L 



P(x, E, L) = x 6M - E + — . (7.2) 

x 6 

Duality acts on M and G as before, and transforms L as 

l ^=(mTTF- < 7 ' 3 > 

The relation ( |5.7| ), apart from the appearance of L and -L as arguments of y and y 
respectively, is unchanged. The earlier treatment can be generalised by defining 

y k = y k (x, E, g, L) = W fc y(w -fc x, ui~ 3Mk E, g, w " 3 ( A/+1 ) fc L) , (7.4) 

so that 



1 d 1 
y x 2 dx x 
If we also define z klk2 as hi ( p. 18 ), then, for k\ and &2 differing by an integer 



Vk - Gi^l s)Vk + e- 2k ™P(x,E,L)y k = . (7.5) 



4"fc 2 - G{~ - \)z klk2 - e- 2k ™P(x,E,L)z klk2 = , (7.6) 

x ax x A 

(with e ~ 2kim = e~ 2k2nl = e~ 2km ) which is the adjoint to ( f7~5|) . Again we can re- 
cover the original problem by shifting k by a half-integer. As before, we find that 
z-i/2,i/2(x, E,g,L) =i\^3y(x,E,g,L) and (cf. ( J2.28Q ) 

T(E, g, L)y_ 1/4 y 1/4 = y_i/ 4 ys/4 + ^-3/4^3/4 + 2/-5/4Z/1/4 • (7-7) 

The presence of a non-vanishing L has however introduced an extra complication. To 
see this explicitly, shift k by ±1/4 to get 



T yoy±i/2 = yoy±3/2 + yw-1/2 + y-m/2 , ( 7 - 8 ) 



_3A/ 



with T = T(Euf^ 4 ,g,=piL). If this equation is rewritten in terms of the function 
y(x,E,g,L), both signs of L appear: for k integer or half-integer we have 

y k = u; k y(xu J - k ,Eu J - 3Mk ,g,(-l) 2k L) . (7.9) 

Notice that the same does not happen for the argument g (or G), which is why this prob- 
lem did not arise before. Prom an analytic point of view, y(x, E, g, L) and y(x, E, g, —L) 
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are just two distinct points of the same function, but in the derivation of the nonlinear 
integral equation it is only the analyticity in E that is used. To proceed, it is best 
to consider L to be held fixed once and for all, and to treat the pair of functions 
Vk = uj k y(io~ k x,io~ 3Mk E,g,L) and v^ = uj k y(uj~ k x,Lj~ 3Mk E, g, — L) independently. 
Then ( |7.8| ) becomes 

r ± U U ± l/ 2 = V0«±3/2 + UlU-l/2 + U-lUl/2 • (7.10) 

This equation is very reminiscent of those given in H for the ^-lattice model. There 
remains an x-dependence in ( [7.10| ) which can be eliminated, once again, by expanding 

v = D + (E,g,L)x+ + D°(E,g,L)xo + D-(E,g,L)x- , (7.11) 

v = D + (E,g,L)x + + D (E,g,L)xo + D-(E ig ,L)x- , (7.12) 

where {x+,Xo>X-} an d {x+,Xo>X-} are alternative bases defined via the behaviour 
near the origin 

Xi(x,E,g,L)~ x x > + 0(x x > +3 ) , i = +,0,-, (7.13) 

X i(x,E,g,L)~x~ x >+0(x~ x > +3 ) , » = +,0,-, (7.14) 

and the Aj's and Aj' are respectively solutions of the indicial equations 

(A - 1)(A(A - 2) - g(g + 2)) + L = , (A - 1)(A(A - 2) - g(g + 2)) - L = . (7.15) 

If the labeling is chosen consistently with that of section §||, so that the A's and the A's 
reduce to the quantities in ( |5.16 ) when L = 0, then 



J- Vo ^±1/2 _ ^0 ^±3/2 ^ ^1 ^-1/2 ' ^-1^1/2 ' v'-lt>) 

with 

Q ± (E,g,L) = E h JM L D ± (E,g,L) , Q± = Q ± {u;- 3Mk E,g, L) , (7.17) 



A±-l 

This leads to two coupled sets of BA equations 



Q ± (E,g,L)=E^M-D ± (E,g,L) , Q± = Q ± {u^ Mk E,g, L) . (7.18) 



°° t?± ,.-3Mr± _ °° p± , ,-3M/2p± 

n * =» - ^A+-A+-i tt h k -" 'k (719) 

00 S± , ,-3Af p± °° z?± , .-3M/2 S± 

Generalising the analysis of §^, it should be possible to derive a nonlinear integral 
equation relevant to this more general case. We expect that this equation will coincide 



with the a 2 -related case of the equations found in [41, 0], in its massless limit, but we 
will leave a detailed investigation for future work. 
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8 Conclusions 

We have continued to study the relationship between integrable quantum field theories 
and ordinary differential equations, and in the process have obtained a novel nonlinear 
integral equation which is able to describe the 0i2, 4>2i and </>is perturbations of minimal 
models within a unified framework. We have also found a natural generalisation of 
the duality symmetry enjoyed by the Schrodinger/massless sine Gordon system [g]. A 
major theme has been that the A2 structures hidden inside certain third-order ordinary 
differential equations, and also inside certain integrable quantum field theories and BA 
systems, are very closely related. It seems clear that the correct way to generalise to yet 
further models is to look to differential equations of even higher order. While this might 
appear to be a task of ever-increasing complexity, there are some reasons to suppose 
that a more unified picture will ultimately emerge. ADE structures have been observed 
in many different, but related, settings in the context of integrable models (see, for 
example, [p|, |23|, |4l"]-45| and references therein). One might hope that the process of 



generalisation will reveal similar phenomena on the differential equations side of the 
correspondence, but more case-by-case analysis will certainly be required before this 
can be confirmed. 
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Notes added — 

(i) The 'massless' nonlinear integral equation derived in §3 has appeared previously, in 



connection with the Izergin-Korepin model, in [46|. 



(ii) A conjecture due to Kausch et al p7, 47| states that the (^12 perturbation of M. Pt 



■■</ 



and the </>i5 perturbation of M. v > ^ have identical ground-state scaling functions if (and 
only if) p' = p/2, q' = 2q. (This implies p = 2 mod 4, since (p, q) and (p',q') must both 
be coprime; such pairs are called 'type II' in p7| .) It is easily checked that this equality 
follows from the recipe for finding ground-state scaling functions given in §6 above: the 
values of M and a that should be used in the two cases are identical, and so both are 
described by the same nonlinear integral equation. We take this as additional support 
both for our conjectures and for that of [37]. 



(iii) In a recent paper [48 1 , Suzuki has independently remarked the relevance of higher- 
order ordinary differential equations to integrable models associated with the algebra 
A n , though with a slightly different emphasis from that adopted above. 
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